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Abstract 

We consider the hadroproduction of Xci and Xc2 mesons and their subsequent 
radiative decays to J /ip mesons and photons in the factorization formahsm of non- 
relativistic quantum chromodynamics, and study the decay angular distributions, 
by means of hehcity density matrices, in view of their sensitivity to color-octet 
processes. We present numerical results appropriate for the Fermilab Tevatron. 
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1 Introduction 

Since the discovery of the J /ip meson in 1974, charmonium has provided a useful labora- 
tory for quantitative tests of quantum chromodynamics (QCD) and, in particular, of the 
interplay of perturbative and nonperturbative phenomena. The factorization formalism 
of nonrelativistic QCD (NRQCD) [1] provides a rigorous theoretical framework for the 
description of heavy- quarkonium production and decay. This formalism implies a separa- 
tion of short-distance coefficients, which can be calculated perturbatively as expansions in 
the strong-coupling constant a^, from long-distance matrix elements (MEs), which must 
be extracted from experiment. The relative importance of the latter can be estimated 
by means of velocity scaling rules, i.e., the MEs are predicted to scale with a definite 
power of the heavy-quark {Q) velocity v in the limit v <^ 1. In this way, the theoretical 
predictions are organized as double expansions in and v. A crucial feature of this for- 
malism is that it takes into account the complete structure of the QQ Fock space, which 
is spanned by the states n — ^s+i^^O -with definite spin 5", orbital angular momentum L, 
total angular momentum J, and color multiplicity C = 1, 8. In particular, this formalism 
predicts the existence of color-octet (CO) processes in nature. This means that QQ pairs 
are produced at short distances in CO states and subsequently evolve into physical, color- 
singlet (CS) quarkonia by the nonperturbative emission of soft gluons. In the limit w — > 0, 
the traditional CS model (CSM) [2] is recovered. The greatest triumph of this formahsm 
was that it was able to correctly describe [3] the cross section of inclusive charmonium 
hadroproduction measured in pp collisions at the Fermilab Tevatron [4] , which had turned 
out to be more than one order of magnitude in excess of the theoretical prediction based 
on the CSM. 

Apart from this phenomenological drawback, the CSM also suffers from severe concep- 
tual problems indicating that it is incomplete. These include the presence of logarithmic 
infrared divergences in the 0{as) corrections to P-wave decays to light hadrons and in 
the relativistic corrections to S'-wave annihilation [5], and the lack of a general argument 
for its validity in higher orders of perturbation theory. While the /c-r-factorization [6] and 
hard-comover-scattering [7] approaches manage to bring the CSM prediction much closer 
to the Tevatron data, they do not cure the conceptual defects of the CSM. The color 
evaporation model [8], which is intuitive and useful for qualitative studies, also leads to 
a significantly better description of the Tevatron data, but it is not meant to represent 
a rigorous framework for perturbation theory. In this sense, a coequal alternative to the 
NRQCD factorization formalism is presently not available. 

In order to convincingly establish the phenomenological significance of the CO pro- 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and vN deep-inelastic 
scattering (DIS), e+e~ annihilation, 77 collisions, and 6-hadron decays may be found in 
the hterature; see Ref. [9] and references cited therein. Furthermore, the polarization of 
■0' mesons produced directly and of J/%1) mesons produced promptly, i.e., either directly 
or via the feed-down from heavier charmonia, which also provides a sensitive probe of 
CO processes, was investigated [10-13]. Until very recently, none of these studies was 
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able to prove or disprove the NRQCD factorization hypothesis. However, HI data of 
ep eJ/ip + X in DIS at HERA [14] and DELPHI data of 77 ^ J/ip + X at LEP2 [15] 
provide first independent evidence for it [16,17]. 

The cross section ratio of Xci and Xc2 inclusive hadroproduction was measured by the 
CDF Collaboration in Run 1 at the Tevatron [18]. The XcJ mesons were detected through 
their radiative decays to J/ip mesons and photons. The J/ip mesons were identified via 
their decay to //"""/x" pairs, and the photons were reconstructed through their conversion 
to e+e" pairs, which provide an excellent energy resolution for the primary photons. The 
increased data sample to be collected during Run 2, which has just started, will allow for 
more detailed investigations of these cross sections, including the angular distributions of 
the XcJ decay products. The latter carry all the information on the XcJ polarization and 
thus provide a handle on the CO processes. 

In this paper, we study the angular distributions of the processes pp Xcj + X 
( J = 1,2) with subsequent decays XcJ ~^ J/ '4'! with regard to their power to distinguish 
between NRQCD and the CSM. Specifically, we consider the polar and azimuthal angles, 
9 and 0, of the J ftp meson in the XcJ rest frame. Once a suitable coordinate system is 
defined in that frame, the full angular information is encoded in the helicity density matrix 
p'ly of the primary Xcj production process, where A and A' denote the Xcj helicities. The 
matrix p'ly, therefore, provides a systematic way of analyzing the complicated details 
of the 6 and dependences. In the following, we study p^ly in four commonly used 
polarization frames as functions of the XcJ transverse momentum pt so as to identify 
optimal observables to discriminate NRQCD from the CSM. 

This paper is organized as follows. In Sec. 2, we list the formulas necessary to evaluate 
the helicity density matrix of pp — > XcJ+X J/^p'j+X for J = 1, 2. In Sec. 3, we present 
our numerical results and discuss their phenomenological implications. Our conclusions 
are summarized in Sec. 4. 



2 Analytic results 

We consider the production processes pp Xcjj + X, where J = 1,2 and j denotes a 
hadron jet, followed by the radiative decays Xcj J li'l in the narrow-width approxima- 
tion. (The presence of the hadron jet, which does not have to be observed, allows for the 
XcJ meson to have finite transverse momentum.) Let \/S be the center-of-mass energy 
of the hadronic collision, M the mass of the XcJ meson, p^ the transverse momentum 
common to the Xcj meson and the jet, y and Uj the rapidities of the latter, A = — J, .... J 
the helicity of the Xcj meson in some polarization frame, and 9 and (j) the polar and 
azimuthal angles of the J/ ip meson in the respective coordinate system defined in the XcJ 
rest frame. Then, the differential cross section can be written as 

E iJl^'^i'CM). (1) 
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where pj^y ^aa' ^^^^^ production and decay processes, respectively, Bj denotes 

the branching fraction of Xcj J/^'li ^'^^ d'^^ = d cos 6 d(f). Invoking the factorization 
theorems of the QCD parton model and NRQCD, we have 

- Xafa/p{xa,l^s)xbfb/p{xb,iif){0^'''[n])— — px\' {ah ^ cc[n]d) . (2) 



dp'rdydy, — — ^ ^ M67r. 

Here, it is summed over the active partons a, 6, d = g,q,q and the cc Fock states n, 
fa/h{xa, Pf) are the PDFs of the beam hadron h, Xa is the fraction of longitudinal mo- 
mentum that a receives from h, jifis the factorization scale, {0^''-'[n\) are the MEs of the 
XcJ meson, and 

Paa' {oh — >• cc[n](i) = ^ T* {ah — > cc[n, \']d) T {ah — >• cc[n, X]d) , (3) 

where it is averaged (summed) over the spin and color states of a and b {d), and T denotes 
the transition-matrix element. We have 

rriT exp(±y) + pr exp(±7/j) 
Xa,b - , l^j 



where uit — ^Jm^^-^ is the transverse mass of the XcJ meson. The partonic Mandelstam 
variables s — (Pa + Pb)^) t — {Pa ~ P)^, and u — (p^ — P)^, where P is the four-momentum 
of the XcJ meson, can be expressed in terms of px, y, and yj as 

s = Pt + + '^PtI^^t cosh(y — 
t = - ptttit cxp(yj - y), 

u = -p^ - pTrriT exp(?/ - y^), (5) 

respectively. Notice that s-\-t-\-u — and sp^ — tu. The kinematically allowed ranges 
of Pt, y, and yj are 

<Pt < 



|?/| < arcosh 



2v^ ' 



y/S -rriT exp(-y) y/S - uit ex.p{y) 

— In < 2/7 < In . (6) 

Pt Pt 

Integrating Eq. (1) over yj, wc obtain the inclusive cross section of pp Xcj + 
Choosing a suitable coordinate system in the XcJ rest frame and defining 



(Ai, A2, e, 0|T| J, A) = T{xcj{X) ^ J/^(Ai, e, 0)7(A2, tt - ^, + tt)). 



(7) 



where Ai = 0, ±1 and A2 = ±1 are the hehcities of the J/tl) and 7 bosons, respectively, 
the decay angular distribution is encoded in the matrix 

u Ea„a.(Ai, A2, e, 0|T| J, Ar(Ai, A2, 9, 0|r| J, A) 

•^Ai,A2 I A1A2 I 

This expression may be conveniently evaluated by observing that [19] 



(Ai,A2,^,0|T|J,A) = f-^DiX_^^{-4>,eA)TU^ (9) 

where Dl^,^{a, (3, 7) = (j, m'\D{a, (3, 7) \j, m) is the representation of the rotation operator 
D{a, (3, 7) = exp(— ^772) exp{—i(3Jy) exp{—iaJz), with a, (3, and 7 being the Euler angles, 
in the eigenstates |j, m) of and Jj. We have 

Din'mi(^,f3,'y) = exp(-z7m')rf^/„(/3) exp(-iam), (10) 

where c?m'm(/^) = (Jj^^'I 6xp(— i/3Jy)|j, m) are the well-known 0? functions, which may be 
evaluated from Wigner's formula 



mm(i+m^i-m') /( j + ^) I ( j _ ^) I ( j + ^/) I _ ^/) I 
f/?! — / 1 \k-m+m' V v-^ ' / w / \J ' / \J / 

"^'"^^^ " .^.axir™-.') ^ - ^ + +m- k)\{j -m'-k)\ 

(^\ 2j+m—m'—2k / ^\ 2k—m+m' 
cosf) (.in^j . (11) 

Owing to the orthogonality relation 

/Att 
d'nDl,^„{-<j>, 9, cP)Di:^.i-<j>, 9, <P) = ^—jSm'm, (12) 



A'ly is normalized as 



J d'nAiy{9,ct>)^5xx', (13) 

so that, upon integration over the solid angle, Eq. (1) reduces to the unpolarized narrow- 
width approximation formula 

^^P^ {PP - Xc.(- J/^7)i + X)^ B{x.j - ^/^7)^^p^ {PP - Xcjj + X) . 

(14) 

By definition, the helicity density matrix pj^^ is hermitian, pft^^ — p'ly. Furthermore, 
invariance of the T matrix under reflection in the production plane, by action of the op- 
erator y = exp{—i7r Jy)V , where V is the parity operator, leads to the symmetry property 
p'L\ _A' = (—1)^"^ Pxx' [19]- Thanks to these constraints, the number of independent real 
quantities contained in p^y is 5 for J = 1 and 13 for J = 2. We may select the entries 
Pxy with A = 0, . . . , J and A' = —A, . . . , A, noticing that they are real if |A'| = A. 
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As for the decay amplitude T'/^;^2, angular- momentum conservation imposes the selec- 
tion rule |Ai — A2I < J. Furthermore, parity conservation entails the symmetry property 
[19] 



[15) 



where rj, rji, and 772 (J, Ji, and J2) are the parity (total-angular-momentum) quantum 
numbers of the XcJ, J/ip, and 7 bosons, respectively. An independent set of T/^a2 ampli- 
tudes thus reads 



tl 


— ^ 1,1 




- ^-1,-1, 


tl 


_7-i 
~ 0,- 


1 


- -T^ 

~ -^0,1' 


tl 


= Tl, 




= T'1,-1, 


tl 


= Tt 


1 


~ -^0,1' 


tl 


rp2 
- ^1,- 


1 


rp2 
- -^-1,1- 



(16) 

Using these ingredients, we can now work out the general form of the angular distri- 
butions of the J/ip meson. 



\,X'=-J 



(17) 



We find 



3 

An 



Po,o 



rQ cos^ 6 -\ — - sin^ 1 
2 



rl sin^^ + ^(l + cos^^) 



^ fplo 
47r 1 2 



^ (1 - 3 cos' 0f + 3rl sin' cos' 9 + ^rj sin^ 9 



+ p1i 



+ P2,2 



2 2 

3r^ sin' 9 cos' ^ + ^(1 - 3 cos' ^ + 4 cos^ 9) + ^ sin' ^(1 + cos' 9) 

\rl sin^ 9^% sin' ^(1 + cos' ^) + ^(1 + 6 cos' 9 + cos^ 9) 
4 2 8 



+ v^Rep^o 
- Repoi 



rl{l - 3 cos' 9) - r\{\ - 2 cos' 9)^^ sin' ^ 



3r^ sin' ^ + 2r\ cos' ^ - ^(3 + cos' ^) 



3r^ cos' ^ + ^(1 - 4 cos') - ^ sin' ^ 



sin ^ cos ^ cos 
sin ^ cos 9 cos 
sin' 6* cos(20) 
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— \/Q Re p\ q 



2 2 

^(1-3 cos^ 9) + rl cos^ ^ - ^(1 + cos^ i 



sin^ 9 cos (20) 



+ Rep|_i \^rl - 2rl + sin=^^cos^cos(3(/)) 
+ %^(^ro^-r? + f)sin^^cos(40)|, 



(18) 



where 







2 




=0 





(19) 



are positive numbers satisfying Z)^=o^a = 1 to be determined experimentally. In fact, 
the Fermilab E835 Collaboration [21] studied the angular distributions of the exclusive 
reactions pp XcJ J li'l ~^ e+e~7, with J = 1,2, and measured the fractional 
amphtudes of the electric dipole (El), magnetic quadrupole (M2), and electric octupole 
(E3) transitions. They found the El transition to dominate, as expected from theoretical 
considerations. From their results, we extract the values 



= 0.498 ± 0.032, 



0.075 ± 0.029, 



0.502 ±0.032, 
0.250 ±0.048, 



0.674 ±0.052. (20) 
= 0.5, rl = 0.1, rl = 0.3, 



The corresponding values for a pure El transition read Tq 
and rl = 0.6. 

We work in the fixed-flavor-number scheme, i.e., we have = 3 active quark flavors 
q — u,d,s va. the proton and antiproton. As required by parton-model kinematics, we 
treat the quarks q as massless. The charm quark c and antiquark c only appear in the 
final state. We are thus led to consider the following partonic subprocesses: 



gg 
gQ 

qq 



cc[n\g, 
cc[n]q, 
cc[n\g. 



(21) 



As for the Xcj mesons, the cc Fock states contributing at LO in v are n — ^Pj^\^Si^\ 
Their MEs satisfy the multiplicity relations 



QXcJ 
QXcJ 



(2J±1)(C»^' 



3d(1) 



3p(l) 

'^Sf^ ) = (2J±l)(C^-« '-^Sl 



3c(8) 



(22) 

which follow to LO in v from heavy-quark spin symmetry. 

Depending on the value of J, the partonic helicity density matrices defined in Eq. (3) 
may be decomposed as 



Paa' 



(ah 



cc 



25+1 r(C) 



4^^] ci)=e;(A)e.(A')F 



cig''' ± C2pX + cM + ^ {p'aPl + pM 
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+ (pX + PlPl) + (pX + PlPa) 



(23) 



where e'^(A) (e^'^(A)) is the polarization vector (tensor) of a J = 1 (J = 2) boson with 
mass M, four-momentum P, and hehcity A, and F and q are functions of the partonic 
mandelstam variables s, t, and u. Since we work at the tree level, F and q are real, so 
that the imaginary parts of pj^^^, all arise from e^(A) and e^^{\). The functions F and 
Ci for processes (21) with n = '^Pi'^M^^ arc listed in the Appendix of Ref. [13]. The 
resuhs for n = ^^f^ may be found in Eqs. (B27)-(B31), (B39)-(B43), and (B51)-(B55) 
of Ref. [11], where one has to identify Ci = a, C2 = 

Lorentz-covariant expressions for e^(A) in four commonly used polarization frames, 
namely, the recoil, Gottfried- Jackson, target, and CoUins-Soper frames, may be found 
in Ref. [22]. These frames differ in the way e^(0) is fixed. In the Xcj rest frame, 
where e^(0) = (0, e(0)), we have e(0) = —p, where p = Pp + p^, in the recoil frame; 
e(0) = Pp in the Gottfried- Jackson frame; e(0) = —pp in the target frame; and e(0) ~ 



{Pp ~ Pp) /y 2 (1 — Pp - Pp) in the CoUins-Soper frame. Here, p = p/\p\ denotes the unit 
three- vector. With the help of the addition theorem for two angular momenta, €^^{X) 
can be constructed from the polarization four- vectors of two J — 1 bosons with the same 
four-momentum as [23] 



where (l,Ai; 1, A2|2,A) are Clebsch-Gordon coefficients. We have P^e''(A) — P^e'^'^(A) = 
gi^ue^^W = and e'^'^(A) = e^'^(A). 

3 Numerical results 

We are now in a position to present our numerical results. We first describe our theoretical 
input and the kinematic conditions. We use rric — M/2 — (1.5 ± 0.1) GeV and the LO 
formula for a^J^^\iir) with nj = 3 [24]. As for the proton PDFs, we employ the LO 
set by Martin, Roberts, Stirling, and Thorne (MRST98LO) [25], with asymptotic scale 
parameter A'^''^ = 174 MeV, as our default and the LO set by the CTEQ Collaboration 
(CTEQ5L) [26], with A^^^ = 192 MeV, for comparison. The corresponding values of A^^) 
are 204 MeV and 224 MeV, respectively. We choose the renormahzation and factorization 
scales to be Hi = ^imT, with i = r, f, respectively, and independently vary the scale 
parameters and ^/ between 1/2 and 2 about the default value 1. We adopt the values 




Y: (l,Ai;l,A2|2,A)e'^(Ai)e^(A2), 



(24) 



Ai,A2=-l 
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of 



3p(l) 



and ( 



^S[^^ ^ from Ref. [12]. Specifically, the former was extracted 
from the measured partial decay widths of Xc2 — ^ 77 [24] , while the latter was fitted to the 
Pt distribution of Xcj inclusive hadroproduction [4] and the cross-section ratio c^xciI'^xa 
[18] measured at the Tevatron. For set MRST98LO, these values read (O^^^ ^P^^ ) — 
(8.9 ±1.3) X 10-2 GeV^ and (O^-'' f^f^]) = (2.3 ±0.3) x 10-^ GeV^. The corresponding 
values for set CTEQ5L are (9.1±1.3) x 10"^ GcV^ and (1.9±0.2) x lO'^ GeV^, respectively 
In order to estimate the theoretical uncertainties in our predictions, we vary the unphysical 
parameters and as indicated above, take into account the experimental errors on rric 
and the default MEs, and switch from our default PDF set to the CTEQ5L one, properly 
adjusting A^^^ and the MEs. We then combine the individual shifts in quadrature, allowing 
for the upper and lower half-errors to be different. 

Our numerical results are presented in Figs. 1-14. Figure 1 is devoted to the cross 
sections a oi pp ^ Xci ± ^ (upper frame) and pp Xc2 + ^ (lower frame). Figs. 2-5 
to the helicity matrix elements p'ly of the former process, and Figs. 6-14 to those of the 
latter one. The matrices pfy are normalized so that their traces, Z]a=-i Paa = "^Pii + Poo 
and J2'x=~2 Px\ = '^P22 ± 2pii ± PoO) unity. We only display the real parts of piy, which 
enters Eq. (18). In each figure, the NRQCD (solid lines) and GSM (dashed lines) results 
are displayed as functions of pt', the central lines indicate the default predictions, and 
the shaded bands the theoretical uncertainties. In Figs. 2-14, four different polarization 
frames are considered: the recoil, Gottfried- Jackson, target, and Collins-Soper frames. 
The results in the Gottfried- Jackson and target frames almost coincide if |A — A'| is even. 
If |A — A'l is odd, the same is true, apart from a relative minus sign. 

We first discuss da/dpr- From Fig. 1, we observe that the GSM contributions essen- 
tially exhaust the NRQCD results at small values of pr, while they get rapidly suppressed 
as the value of pr increases. This may be understood by observing that, with increasing 
value of Pt, the processes (21) with n = '^S[^^ gain relative importance, since their cross 
sections involve a gluon propagator with small virtuality, = M^, and are, therefore, 
enhanced by powers of p^/M'^ relative to those of the other contributing processes. In the 
fragmentation picture [27], these cross sections would be evaluated by convoluting those 
of gg — > gg, gq — > gq, and qq gg with the g ^ cc ^S^^ fragmentation function [28]. 
For such processes, the attribute fragmentation prone has been coined [29]. 

We now turn to dp\y/dpT. Looking at Figs. 2-5, we observe that the NRQGD and 
GSM predictions are generally rather similar in the low-pT range. However, at large values 
of Pt, there may be dramatic differences, depending on the matrix element p\y considered 
and the polarization frame chosen. Specifically, the diagonal elements, Pq q and p]^^, lend 
themselves as powerful discriminators between NRQCD and GSM in all four polarization 
frames. In the case of p\q, only the Gottfried- Jackson and target frames are useful, while 
the GoUins-Soper frame is preferable in connection with p\ 

We finally draw our attention to dp^y/dpr (see Figs. 6-14). Again, the NRQCD and 
GSM predictions always merge in the limit px 0. On the other hand, in the large- 
Pt regime, we can always find a polarization frame that allows us to clearly distinguish 
NRQCD from the GSM. As for pi and pli, all four frames can serve this purpose. 
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The Gottfried- Jackson and target frames work best for pIq, pl^^, and pl_i, while the 
Collins-Soper frame is the frame of choice for Pi_i, pi 2) P2 0> P2-2- We observe 
that the NRQCD results for pf^, P20) P2-I) ^''^^ P2,-2 numerically suppressed, with 
magnitudes of order 0.1 or below. Notice that the NRQCD and CSM results for pl^, 
(A' = 2, 1, 0, —1, —2) differ, although the only contributing cc Fock state with J = 2 is 
a CS state, n — ^P2'^. This may be understood by observing that the CO contribution 
enters through the normalization of p\y. 



4 Conclusions 

Run 2 at the Tevatron, which has just begun, will provide us with a wealth of new 
information on how cc pairs turn into physical charmonia. In particular, this will allow 
us to put the NRQCD factorization hypothesis to a stringent test. In the face of this 
exciting situation, we considered the processes pp — > XcJ + ^ for J = 1,2 followed by 
the radiative decays XcJ ~^ J/ '4'! under Tevatron kinematic conditions in the NRQCD 
factorization formalism, and we investigated the decay angular distributions, using the 
helicity density matrix formalism, with regard to their power to identify the presence of 
CO processes. Specifically, we expressed the distributions in the J/^jj polar and azimuthal 
angles in the XcJ rest frame in terms of the helicity density matrix elements pxy the XcJ 
production processes. We then analyzed the matrix elements p'ly as functions of the XcJ 
transverse momentum in four frequently used polarization frames, namely, the recoil, 
Gottfried- Jackson, target, and Collins-Soper frames. We found that the CS processes play 
the leading role in the low-pr range. This fact could be exploited to extract the CS ME 



QXcO 



3p(l) 



At large values of pr, typically for pt ^ 5 GeV, the NRQCD and CSM 

predictions can differ significantly, depending on the matrix element p'ly considered and 
the polarization frame selected. We could demonstrate that for every matrix element p'ly 
there is at least one polarization frame where the NRQCD and CSM results are greatly 
different. The decay angular distributions under consideration here should, therefore, 
lend themselves as sensitive probes of the CO processes to be doubtlessly estabhshed by 
experiment. 
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Figure 1: Differential cross sections da/dpT of pp Xcj + ^ for J = 1 (upper frame) 
and J = 2 (lower frame) at = 1.8 TeV in NRQCD (solid lines) and the CSM (dashed 
lines). The central lines represent the default predictions, and the shaded bands indicate 
the theoretical uncertainties. 
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Figure 2: Differential helicity density matrix element dpl ^/dpT of pp — > Xcj + X at 
\/S = 1.8 GeV in four different polarization frames in NRQCD (solid lines) and the GSM 
(dashed lines). The central lines represent the default predictions, and the shaded bands 
indicate the theoretical uncertainties. 
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Figure 7: Same as in Fig. 2, but for dpi i/dpT- 
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Figure 8: Same as in Fig. 2, but for dRe pl^/dpT- 
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Figure 9: Same as in Fig. 2, but for dpi _i/dpT- 
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Figure 10: Same as in Fig. 2, but for dpl2/dpT- 
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Figure 11: Same as in Fig. 2, but for dUe P21I dpT- 
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Figure 12: Same as in Fig. 2, but for dUe p^^l dpT- 
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Figure 13: Same as in Fig. 2, but for dRepl _i/dpT- 
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Figure 14: Same as in Fig. 2, but for dp^ _2/(ipT- 
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